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O ' 

We study a recent generalization proposed for the XY model in two and three dimensions. Using 
\ both, the continuum limit and discrete lattice, we obtained the vortex configuration and shown 

that out-of-plane vortex solutions are deeply jeopardized whenever the parameter of generalization, 
L, is increased. The critical temperature for such models is calculated using the self consistent 
harmonic approximation. In both, two- and three-dimensional cases, such a temperature decreases 
with raising L. Our results are also compared with other approximated methods available in the 
• literature. 
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INTRODUCTION AND MOTIVATION 



The Hamiltonian H = —JJ2<ij>(S? $j + $f Sj), where J is a coupling constant, the summation is taken over the 
nearest-neighbors sites in a square lattice and S" (a = x, y) is a spin component at site i, describes two important 
models in magnetism known as the two-dimensional (2D) planar rotator model (PRM) and the XY model. The 
difference between these two models lies in the number of spin components. Thus, in the PRM framework we have 
only two spin components which are subject to the constraint S 2 +S 2 = S 2 . On the other hand, XY model also displays 
the third component with the condition S 2 + S 2 + S 2 = S 2 . Hence, whenever fixing the spin 5* of the system, we are led 
with one and two independent spin degrees of freedom for these models, respectively. Thus we may parametrize the 
' PRM physical spin content by only one scalar field, the azimuthal angle ipt, say, (Si)pRM = S(cosipi, siinpi), while in 
the XY model, we have two scalar, the azimuthal and polar angles, say, (Si)xY — iS(sin#i costpi, sin^ sint^, cos#i). 
As it is well-known, the 2D PRM and XY models support topological excitations and, although no long range order is 
established at any finite temperature, they are shown to exhibit phase transition related to the unbinding of vortex- 
antivortex pairs. Indeed, at low temperatures, the vortex excitations are bound into pairs whose interaction grows 
logarithmically with the distance between the vortex centers. However, whenever the so-called Berezinskii-Kosterlitz- 
' Thouless (BKT) transition temperature, Tbkt, is reached, the bound pairs appear to dissociate. 

In addition, the static critical behavior of such models, mainly the PRM, in two and three dimensions, has been 
studied for a long time, using several statistical and Monte Carlo (MC) simulation methods. Among others, such 
studies have led to the consensus upon the nature of the phase transitions and the values of the critical temperature 
and exponents^. On the other hand, an interesting extension of the 2D PRM was introduced around two decades ago 
by Domany, Schick, and Swenden^, in which the potential is given by V(x) — 2.7(1— [cos 2 (x/2)] p ), where \ = <pj — <pi. 
Clearly, the usual PRM is recovered as long as p = 1. Nevertheless, as p is raised, an increasingly narrow potential 
well, with a width around n/p, is observed in the model. Notice also that whenever we have x — (^n + 1)tt (n 
■ integer), then the potential gets a constant value, V((2n + l)~ir) = 2 J. In addition, ordinary spin-wave behavior at low 
temperatures is ensured by the fact that, near the minimum we do have V(x) ~ Jp 2 X 2 /2- It was also shown by MC 
simulation that, for very large p a first-order, rather than a continuous phase transition takes placed. In this sense, 
it seems that the vortex excitations alone are sufficient to account for both continuous and first-order transitions, 
H though in qualitatively distinct wayai. 



Recently, a generalization of XY model has been proposed by Romano and Zagrebnov 5 , whose potential looks 
like — J(sin6*i sin6j) L cos(ifi — ipj) , with L G N. In their paper— the authors have established rigorous inequalities 
holding for every L > 1. Using Mean Field (MF) and Two-Site Cluster (TSC) techniques they have also estimated the 
critical temperature, in the 3D case, for some values of L. In the 2D case, and for arbitrary L, they suggested that the 
above potential produces orientational disorder at all finite temperatures, and supports BKT transition. However, the 
stability and behavior of the vortex-like solutions in these models are still unclear. Therefore, it would be important 
to know about the types of vortices present in the system and their roles in such a phase transition. Thus, in this work 
we present a detailed analysis of the stability of vortex solutions. This is done by applying discrete lattice and usual 
continuous approximation approaches. The main conclusion is that out-of-plane solutions are deeply jeopardized as 
long as the anisotropic parameter, L, is increased. Furthermore, concerning the problem of phase transitions it would 
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also be desirable to verify in more details the validity of some approximated methods used to estimate the values of 
the critical temperatures. For this purpose, we apply here the Self Consistent Harmonic Approximation (SCHA) to 
the generalized XY model in two and three dimensions. Our results strongly supports the idea that as L is turned up 
the critical temperature of the 2D and 3D generalized XY models appears to decrease. However, it is worth noticing 
that our results predict lower critical temperatures than those obtained by both, MF and TSC methods. Namely, 
for the case L = 1, SCHA approach yields the value which is the closest to that found by means of Monte Carlo 
simulations. 

The generalized XY model is defined like below: 



,- X) Y, [i-(s?) 2 -(s*) 2 + (s*stf} L/ \s*s* + s!s^ 



TjGen j 
Wyv — — J 

<ij> 

= — J \_. ( sm #i s'm6j) L cos(cfi — ipj) , (1) 



where the spin vector at site i is defined by «Sj = (Sf , Sf, Sf) satisfying the non-linear constraint Sf = 1, and L e N. 
Clearly, as long as L = 1 we recover the usual XY- model, while for L = we get the PRM (with S z — 0). In addition, 
as long as L increases the anisotropy of the generalized model is taken up. This Hamiltonian is a function of the spin 
components of the form f(Sf,S z 3 )[SfS] + SfS]} 7 where f(Sf , S?) = -[1 - (Sf) 2 - (S*) 2 + (SfS]) 2 ] L ' 2 . 



II. THE SOLUTIONS AND THEIR STABILITY 



In order to study possible vortex-like and spin-wave configurations, as well as their stability, we need to take into 
account the continuous version of Hamiltonian It is not difficult but length to obtain that: 



jrGcn 

ti XY - 



J 



(1 - m 2 ) L (W) 2 + L 2 m 2 (l - m 2 ) L - 2 (Vm) 2 - 4(1 - ™ 2 ) L + 4 



(2) 



where to = cos#, while the constant +4/a 2 has been introduced in order to renormalize the energy. 

The equations of motion for spin-field variables, m(x,y;t) — cos[9(x,y;t)] and <p(x,y;t) read like follows: 

J(l - m 2 ) L - 1 |2mZ(Vm) ■ (V0) - (1 - to 2 )V 2 <; 

L(Vm) 2 L(L — 2)m 2 (Vm) 2 LmV 2 m 



dm 
~dt 



Of 



JLm(l — m 



1\L-1 



(1 — rn 2 ) 



(1-to 2 ) 2 



(1 — m 2 ) 



(3) 
(4) 



Now, we shall pay attention to the possible topological excitations associated to the non-linear equations of motion, as 
well as their stability. First of all, we notice that genuine static planar vortex-like configurations, to = cos 9 = and 
4> = arctan(j//x), are supported by the generalized model, for every L (L > 0, L integer). This can be easily checked in 
Eqs. © and Later, let us notice that whenever rrii = 0, then sin#i = +1, implying that L is immaterial, since the 
energy experiences no changes as L is varied. Nevertheless, let us suppose quasi-planar configurations, say, rrii ps and 
sin6*i < 1. Now, the scenario is deeply changed: indeed, increasing L yields to smaller effects of the term (1 — to 2 ) l (or 
f(Sf, Sj)) in the Hamiltonian, in a such a way that, as L becomes very large, this term practically vanishes, leading to 
an appreciable decreasing of the spin interactions, and an increasing of the system energy. Therefore, we may conclude 
that, as long as the parameter L is taken up, the out-of-plane spin components of the generalized XY-model appear 
to become more absent in the system. In other words, such components would require a greater amount of energy 
to show up. It indicates that out-of-plane vortex solutions must be unstable. The fact that out-of-plane fluctuations 
((^i) 2 ) must decrease as L increases will be used to calculate the critical temperature by means of the SCHA. 

Now, in order to analyze the stability of possible out-of-plane vortex solutions in more detail, we use the technique 
which has been applied by Wysin^ to determine the critical anisotropy parameter, A c , for the anisotropic Heisenberg 
model 6,7,8 , whenever out-of-plane vortex solutions become unstable^. This technique is based upon the discrete 
lattice, and consists in analizing the energy of out-of-plane vortex solutions, along with non-zero out-of-plane spin 
components appearing only in a finite number of sites equidistant from the vortex center. In a first approximation, 
we consider that only the first four nearest spins of the vortex center have non-zero out of plane components, mi =/= 0, 
while all the other sites display to = 0. In order to improve this calculation, we can include more sets of equidistant 
sites. Considering the first three sets of equidistant spins at radius r\ = 1/V2, r 2 = VlO/2 and r 3 = \/l8/2, with 
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out-of-plane components toi , to 2 and to 3 respectively, we then find the effects due to the 40 bond sites nearest to the 
vortex center: 



E 3 = -16 J 



(1 - m 2 ) L/2 , 4(1 - m 2 ) L/ ' 2 , (1 - m\ + m x m 2 - to 2 ) l / 2 

(5) 



5 V5 \fl 

^ _ m 2)L/2 + (l~m 2 +m, 2 m 3 ~m 2 ) L / 2 5 



3 V J ' V5 2 

where we have subtracted out the ground state energy —40 J (in- plane exchange). 

Minimization of the energy with respect to the out-of-plane spin components demands the vanishing of the following 
derivatives: 

JL -(m 2 -2mi)(l-m\ + mim 2 -ml) L/2 ^ 1 1 (6a) 



dmi V5 

dE_ _ _^£( mi _ 2 TO2 )(l - m f + mi m 2 - mj) 1 ^^ 1 + ^-^m 2 (l - m 2 ) L/2_1 



dm 2 ^5 - 

il4t ^TO 2 (l - m 2 ) L/2_1 - ^^(m 3 - 2m 2 )(l - to| + to 2 to 3 - to 3 ) l/2 ~\ (6b) 



V65 V5 

9-E 8JL n \/ 2 2\L/2— 1 

= -j^\m2 - 2to 3 )(1 - m 2 + m 2 m 3 - m 3 ) ' 

Om 3 v5 

+ 16JLy|m 3 (l-m 3 ) L/2 - 1 . (6c) 

Then, out-of-plane vortex solutions are stable whenever Eqs.© identically vanish, while rrii ^ 0, for every i = 1,2, 3. 
To solve these equations, we assume small amplitudes for rm. In this case, the linearized version of these equations 
reads like below: 

16JL 8JL 

-yg-wi - "y|- m 2 = °> ( 7a ) 
8JL /16JL 32 JL MJL\ 8JL 

"7T mi + l^ + ^r + wJ m2 "7T m3 = ' (7b) 

- ^m 2 + 16 JL + \/i}m 3 = 0. (7c) 

Of course, trivial solutions, (rrii — 0), are always possible (in-plane vortex). Nevertheless, our interest lies in 
the conditions which could provide non-trivial configurations to show up. The latter ones are possible if and only 
if the determinant of the matrix coefficients is zero. However, such a determinant takes the following value: d ss 
3331. 52( JL) 3 , which implies that non-trivial configurations are possible only for L = (PRM, which of course, can 
not support out-of-plane solutions). This result clearly states us that out-of-plane vortex solutions are unstable. 
Furthermore, we have also performed a similar plain for equidistant sets of spins from r\ = 1/V2 to r 6 = V50/2 
(six sets of spins) , whose conclusion is analogous to the above one. In Ref . 6] , Wysin found an accurate value of A c 
(critical anisotropy) for the anisotropic Heisenberg model which is in good agreement with Monte Carlo results, using 
only three sets of equidistant spins (form r\ = l/\/2 to r 3 = \/l8/2). Thus, we claim that as long as we have used 
six sets of equidistant spins, our results can be generalized leading us to conclude that out-of-plane vortex solutions 
will be unstable in this model, JQ), for every L value. 

Analogous conclusion is also obtained if we study the plane-wave (meson) behavior of the model. In order to see 
that, let us take expressions !• !l II) . and suppose the appearance of small time-dependent deviations from the planar 
vortex-like solution (m = too = cos0 = and (f> = 4>q = arctan(y/x)), like below: 

to — > m = too + toi (x; t) = TOi (x; t) , (8) 
4> — > 4> = 4>o + 4>i (x; i) = arctan(y/a;) + 4>i(x;t) . (9) 
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Since mi and <p± are small, we may neglect higher orders terms proportional to them and their derivatives in the 
equations of motion-. In this case we would be neglecting the spin wave interactions. In addition, taking mi(x; t) = 
mi(x) e luJt and 4>i(x;t) = (f>i(x) e tuJt , where u> is the frequency of such long wavelength spin-wave excitations, we get 
the following equations for mi and 4>\\ 

V 2 <Mx) = -ymi(f), (10) 

iui 

which imply that: 

V 2 0i(x) = - f _ M*) • (12) 

Then, as r — > oo, we expect that such excitations look like plane-waves, say, m\{x) — m\ e tk ' x and </>i(x) = 0" e 4 ^, 
getting: 

+4/ 2 

c 2 (|fc|;L) = -^|fc| 2 L, (13) 

which confirms our earlier conclusion. Hence, for L > 1 out-of-plane fluctuations appear to spend more and more 
energy. As L is increases, the magnon density decreases as e^^, while the planar- vortex density is not affected. It 
may imply in important differences in the thermodynamics of these Generalized XY systems. In fact as has been 
shown by Currie et. al. 10 , the phase shift interaction between spin waves and solitons provide the sharing mechanism 
of energy and degrees of freedom among the nonlinear excitations of the system and therefore is important in the 
study of the statistical mechanics of the model. 



III. PHASE TRANSITIONS IN TWO- AND THREE-DIMENSIONAL GENERALIZED XY MODELS 



Romano and Zagrebnov in Ref. p| have shown that this model supports Berezinskii-Kosterlitz-Thouless (BKT) 
phase transition in 2D, and ordering transitions in 3D. In this section, we perform an estimative of the critical 
temperature, in 2D and 3D, by applying the Self Consistent Harmonic Approximation (SCHA), which has been 
widely appeared in the Iiteratureii*i2ii2ii£ii&ii&, in connection with isotropic and anisotropic models of the same class 
of universality of the XY model. 

Such a technique works by replacing the original Hamiltonian of the model by a harmonic one, along with the 
coupling constant, J, by an effective one, K, which takes into account the nonlinearities of the interactions. First, we 
shall apply this technique to the spatially anisotropic 3D generalized XY model, defined by the following hamiltonian: 



H%n — — J 



^[l-(5 r 2 ) 2 -(5 r \ a ) 2 + (^^ +a ) : 



L/2 



cos((p r 



-J z £ [l - (S*) 2 - (5 r 2 +c ) 2 + {S?;S' i+c f 



L/2 



COs(^ r - <yS r +c), 



(14) 



where r + a and r + c represents the nearest-neighbors in the XY plane, and in the z direction respectively, and J z 
is the inter-plane coupling constant. The harmonic hamiltonian related to the hamiltonian i|14|) is given by: 



r a c 

+ (4LJ + 2LJ Z )(S*) 2 



where K and K z are given by: 

K = j([l-(S r *) 2 -(S* +a ) 2 - 
^ = J*([l-(S*) 2 -(S r * +c ) 2 



(^5 r % a ) 2 ] L/2 cos(^-^ r+a )), 



(^ +c ) 2 ] i/2 cos(^ r 



■)>■ 



(15) 

(16) 
(17) 
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These averages are approximated replacing (. . .) by (. . .)o, i.e., they are calculated using Hq instead H. In this case, 
ip T and are decoupled allowing us to write: 



J((l-(^) 2 ) L )(cos(^-^ r+a )), 
K z » J 2 ((l - (5 r z ) 2 ) i )(cos(^ r - Vr+C )). 



(18) 
(19) 



As we have seen in Section II, as L increases, the out-of-plane spins fluctuations decreases. Using this result, we 
expand (1 — (Sf) 2 ) L as 1 — L(S*) 2 + . . ., keeping only the terms of the order of (S*) 2 . Then, the first average is easily 
solved givingMii^: 



1 - L{{St 



= 1 - 



T 



-L, 



(20) 



4 J + 2 J z 

Hamiltonian (|15fl . except for the L term in (S 2 ), has the same form of the Hamiltonian of Ref. |15|, So, for a generic 
value of L, we obtain, in the two limits J z ~ J and J z <C J, the following expressions: 
(a) J Z ^J 



K = J 1 



T 



4 J + 2 J z 





T 


) expj 


~16 



3if + K Z 



3K + K Z 



1 



2if VI -.9 



2K 2 + KK Z 2K 2 + 2KK 
1 ff tanh _1 (l - gfl 2 



(21) 



where g = K/K z 
(b) J, « J 



ifx = J, 1 - 



T 



4 J + 2J Z 



L ' eXP< > _ 16 



2X + 2JT 2 



2A' + 



if 2 + 2^if z (if + K z f 



1 / tanh x (l - g) 



1/2 



if«J 1 



4 J + 2 J z 



L exp 



K\ (I -gf' 2 



T 
16K 



1-9 



1 



2 l + K z /K 2 + K z /K 



12\ /if 



24tt 



if 



^ In ^ 



if 2 



if 



(22) 



(23) 



if , « J z 1 



T 



4 J + 2J Z 



if 



Td/16X 



exp 



T 
16K 



1 



1 + K Z /K 2 + K z /K 



(24) 



where d= (tt 2 + 12)/12tt and c = (tt/12) ln(4/?r) + (1/tt) ln?r. 

These self-consistent equations for the coupled variables if and if z gives the effective coupling at each value of 
the temperature. The critical temperature is reached whenever these equations admit only the trivial solutions, say, 
if = and if z = 0. In this case, numerical data obtained from these equations can be approximately fitted by 



2\ 



e + L' 



(25) 



where A = J z /J. 

i — in 

For L — and L — 1 we recover previous results from Refs. |15U16| . for the PRM and XY models, respectively, 
which are known to be in good agreement with those obtained via Monte Carlo. In table HI results from Monte Carlo 
simulations for the PRMii and XY models^, and those obtained by Romano and Zagrebnov^ via MF and TSC 
techniques are compared with those calculated by SCHA, for the isotropic case. It is worth noticing that our results 
generally yields critical temperatures lower than those predicted by MF and TSC techniques. Moreover, for the L = 1 
case (usual XY model) for which we have a more detailed study, our result is the closest to that obtained by MC 
simulation. 

On the other hand, whenever J z — 0, we get the 2D XY model. In this case, the effective coupling constant, if, 
representing the stiffness or superfluid density p s , is given by 



(26) 
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TABLE I: Mean-Field, Two-Site Cluster, Monte Carlo and SCHA results for the critical temperature of the isotropic 3D 
generalized XY model. 



L Monte Carlo Mean Field Two-Site Cluster SCHA 

2.2 ±0.05 2.2073 

1 1.54 ±0.01 2.0000 1.7130 1.6136 

2 1.6000 1.4389 1.2716 

3 1.3741 1.2827 1.0493 

4 1.2190 1.1779 0.8931 




FIG. 1: Stiffness as a function of the temperature for L = 1 (a), L = 2 (b), L = 3 (c) and L = 4 (d) in the two dimensional 
case. 



Similarly, the critical temperature is obtained whenever the self-consistent equation for K is solved only by taking 
K = 0. Then, we found that the critical temperature, as function of L, is T = 4/(e ± L), which is quite above 
Monte Carlo results for PRM and XY modelsiS^. This occurs because the SCHA approach attributes an excessive 
energetic cost to topological excitations, such as vortex, causing an overestimation of the Berezinskii-Kostcrlitz- 
Thouless transition temperatur e 11 ! 15 . In Fig. we show the stiffness K as a function of T for some values of the 
parameter L. 
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IV. CONCLUSIONS 



In summary, we have obtained the continuum limit for the Generalized XY- model, defined by Hamiltonian (1) and 
studied the vortex stability. Calculations using the discrete lattice confirmed the results obtained by the continuum 
approach. Out-of-plane fluctuations, and as a consequence, the magnon density, decrease as the parameter of gen- 
eralization L increases. The phase transition temperature as a function of L was estimated using the self consistent 
harmonic approximation. Our results seen to be more realistic than other approximated methods available in the 
literature, such as two-site-cluster and mean field. However, it has to be shown numerically and is beyond the scope 
of this letter. We would like to consider in a future paper some Monte Carlo simulations to calculate the critical 
temperature for the Generalized XY-model in order to confirm our assertion. 
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